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Constant  temperature  molecular  dynamics  has  been  used  to  simulate  the  adsorption  of  hy¬ 
drated  halide  ions  X'  =  F',  Cl',  Br  and  1',  and  lithium  ion  Li*  on  a  flat  uniformly  charged 
surfaces.  The  simulations  were  done  with  either  214  water  molecules  and  two  ions  (l.i+  and 
X')  in  a  box  2.362  nm  deep  or  with  430  water  molecules  and  the  two  ions  in  a  box  4.320 
nm  deep.  The  boxes  were  periodically  replicated  in  the  xy  directions.  The  magnitude  of  the 
surface  charge  on  the  box  ends  was  ±0.11  c/(nm)\  corresponding  to  an  electric  field  of  2 
xlO7  V/cm.  The  lateral  dimensions  of  the  simulation  cell  were  1.862  nm  x  1.862  nm  (x  x 
y)  in  each  ease.  All  of  the  water  molecules  and  ions  interacted  with  the  end  walls  via  a  weak 
9  -  3  potential.  The  ST2  water  model  and  parameters  optimi/.cd  lor  alkali  halides  interacting 
with  the  model  ST2  water  molecule  were  used  in  the  calculations.  Common  practices  of 
truncating  the  interactions  at  a  finite  distance  (0.82  nm)  and  switching  off  Coulomb  inter¬ 
actions  at  small  distances  were  followed.  The  temperature  was  set  at  T  =  2.411  k.I/mole  (290 
K). 

Some  of  the  properties  calculated  were:  distribution  density  profiles  for  ions  and  water  across 
the  gap  important  for  comparisons  with  Gouy-Chnpman  theory,  adsorbed  ion  -  water  pair 
correlation  functions,  the  number  of  water  molecules  in  the  first  and  second  hydration  shells 
of  the  ions  as  a  function  of  time.  The  time  spent  by  a  water  molecule  in  the  hydration  shell 
was  calculated  to  be  approximately  ten  times  longer  for  lithium  than  any  other  ion.  The 
correlation  between  distance  from  the  electrode  and  hydration  number  was  studied  and  gen¬ 
erally  found  to  be  pronounced  for  the  larger  anions.  Comparison  of  the  dynamics  of  the 
common  ion  Li*  for  different  anions  revealed  the  subtle  influence  of  a  transccll  interaction 
in  the  2.362  nm  thick  film. 

In  the  given  field  the  smallest  ions  Li*  and  F'  remained  fully  solvated  at  all  times.  Chloride 
behaved  quite  differently.  Part  of  the  time  this  ion  was  far  enough  away  from  the  electrode 
to  be  fully  hydrated  and  part  of  the  time  it  was  in  physical  contact  (ic,,  phystsorbed)  on  the 
electrode  with  no  water  molecules  interposed  between  it  and  the  electrode.  Bromide  favoured 
contact  adsorption  over  full  hydration  most  of  the  lime.  Iodide  was  observed  to  be  contact 
adsorbed  almost  all  of  the  time.  These  simulations  provide  new  insights  on  the  behavior  of 
strongly  hydrated  ions  at  surfaces  and  how  the  transition  from  non-contact  to  contact' 
adsorption  occurs. 
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ABSTRACT 

Constant  temperature  molecular  dynamics  has  been  used  to  simulate  the  adsorption  of  hy¬ 
drated  halide  ions  X*  =  F*,  Cl*,  Br*  and  1*,  and  lithium  ion  Li+  on  a  flat  uniformly  charged 
surfaces.  The  simulations  were  done  with  either  214  water  molecules  and  two  ions  (Lif  and 
X")  in  a  box  2.362  nm  deep  or  with  430  water  molecules  and  the  two  ions  in  a  box  4.320 
nm  deep.  The  boxes  were  periodically  replicated  in  the  xy  directions.  The  magnitude  of  the 
surface  charge  on  the  box  ends  was  ±0.11  c/(nm)2,  corresponding  to  an  electric  field  of  2 
xlO7  V/cm.  The  lateral  dimensions  of  the  simulation  cell  were  1.862  nm  x  1.862  nm  (x  x 
y)  in  each  case.  All  of  the  water  molecules  and  ions  interacted  with  the  end  walls  via  a  weak 
9-3  potential.  The  ST2  water  model  and  parameters  optimized  for  alkali  halides  interacting 
with  the  model  ST2  water  molecule  were  used  in  the  calculations.  Common  practices  of 
truncating  the  interactions  at  a  finite  distance  (0.82  nm)  and  switching  olT  Coulomb  inter¬ 
actions  at  small  distances  were  followed.  The  temperature  was  set  at  T  =  2.411  kJ/mole  (290 
K). 

Some  of  the  properties  calculated  were:  distribution  density  profiles  for  ions  and  water  across 
the  gap  important  for  comparisons  with  Gouy-Chapman  theory,  adsorbed  ion  -  water  pair 
correlation  functions,  the  number  of  water  molecules  in  the  first  and  second  hydration  shells 
of  the  ions  as  a  function  of  time.  The  time  spent  by  a  water  molecule  in  the  hydration  shell 
was  calculated  to  be  approximately  ten  times  longer  for  lithium  than  any  other  ion.  The 
correlation  between  distance  from  the  electrode  and  hydration  number  was  studied  and  gen¬ 
erally  found  to  be  pronounced  for  the  larger  anions.  Comparison  of  the  dynamics  of  the 
common  ion  Li  +  for  different  anions  revealed  the  subtle  influence  of  a  transcell  interaction 
in  the  2.362  nm  thick  film. 

In  the  given  field  the  smallest  ions  Li  +  and  F  remained  fully  solvated  at  all  times.  Chloride 
behaved  quite  differently.  Part  of  the  time  this  ion  was  far  enough  away  from  the  electrode 
to  be  fully  hydrated  and  part  of  the  time  it  was  in  physical  contact  (ic.,  physisorbed)  on  the 
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electrode  with  no  water  molecules  interposed  between  it  and  the  electrode.  Bromide  favoured 
contact  adsorption  over  full  hydration  most  of  the  time.  Iodide  was  observed  to  be  contact 
adsorbed  almost  all  of  the  time.  These  simulations  provide  new  insights  on  the  bchavioi  of 
strongly  hydrated  ions  at  surfaces  and  how  the  transition  from  non-contact  to  contact 
adsorption  occurs. 
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{.INTRODUCTION 

In  a  previous  paper  (I),  constant  energy  molecular  dynamics  simulations  of  the 
adsorption  lithium  and  iodide  ions  from  aqueous  solution  onto  charged  surfaces  were 
described.  This  paper  continues  this  study  with  an  investigation  of  the  halide  series 
(fluoride,  chloride,  bromide,  and  iodide)  characterized  by  increasing  ionic  radius  and 
decreasing  hydration  enthalpies,  adsorbed  on  an  idealized  flat  charged  surface  situated 
at  the  ends  of  a  box  replicated  periodically  in  the  lateral  (xy)  directions. 

Capacitance,  chronocoulombic,  spectroscopic  and  cllipsomctric  measurements  of  ions 
adsorbed  at  elcctrochmcical  interfaces  have  long  been  interpreted  using  a  standardized 
model  (2-4).  It  has  been  shown  that  the  larger  ions  adsorb  with  at  least  partial  loss 
of  the  inner  solvation  shell.  For  metal  electrodes  this  process  is  called  contact 
adsorption,  and  on  some  noble  metals  a  series  has  been  determined  in  which  higher 
members  displace  lower  ones  from  the  surface  (2,  .5).  A  fuller  discussion  was  given 
in  the  first  paper  (1).  It  is  a  goal  of  this  paper  to  understand  more  about  the  tran¬ 
sition  from  adsorbed  but  fully  solvated  ion  to  adsorbed  ion  in  physical  contact  with 
the  surface  when  there  is  a  static  electric  field  attracting  the  ion  to  the  surface.  We 
tacitly  assume  this  local  property  of  the  ions  is  not  greatly  affected  by  gross 
misshandling  of  the  long  range  electrostatic  fields  of  the  interface  through  the  appli¬ 
cation  of  the  cut-off.  The  field  chosen  is  at  the  upper  end  of  fields  in 
thermodynamically  stable  double  layers.  In  weaker  fields  the  average  position  of  the 
ion  was  closer  to  film  center  (6).  We  show  how  essential  features  of  the  contacting 
process  can  be  modeled  with  a  system  consisting  of  two  ions  and  214  water  molecules 
sandwiched  between  charged  nonmetallic  plates  2.302  nm  apart.  More  elaborate  cal¬ 
culations  with  430  waters,  two  ions  and  plates  4.320  nm  apart  were  performed  to 
check  for  transcell  interactions  between  adsorbed  ions.  Such  interactions  were  indeed 
found.  The  study  revealed  new  features  about  the  way  the  adsorption  process  pro¬ 
ceeds.  Properties  calculated  and  reported  here  include:  time  dependence  of  the  z 


coordinates  of  ions,  time  dependence  of  ion  hydration  numbers,  time  averaged  dis¬ 
tribution  of  ions  and  water  across  the  electrolyte  film,  and  the  ion-water  pair  corre¬ 
lation  functions.  The  first  property  showed  that  the  transition  to  contract  adsorption 
first  occurs  for  Cl  in  the  given  electric  field.  'Hie  second  showed  how  solvation  shells 
were  affected  as  the  adsorbed  ions  jittered  back  and  forth  under  opposing  influences 
of  thermal  forces  and  the  combined  attraction  of  the  electric  and  wall  forces. 

All  the  frequently  used  water  models  (ST2,  TIPS  and  SPCH)  (7-h),  have  parameters 
optimized  to  bulk  properties.  Consequently  none  may  be  adequate  to  describe 
quantitatively  a  water  molecule  outside  a  mean  field  approximation  to  it's  bulk  envi¬ 
ronment.  Our  strategy  to  side  step  this  difficulty  is  to  perform  calculations  on  a  series 
of  related  ions  in  the  expectation  that  by  considering  a  series  with  the  same  charge 
and  a  graduated  property  like  radius,  the  watcr-ion-wall  interactions  will  span  the 
experimental  range  of  phenomena  observed  as  ions  approach  a  charged  surface.  The 
presence  of  the  Li '  counter  ion  in  ail  the  calculations  provides  a  useful  sanity  check 
against  pathological  conditions  that  might  arise  due  to  the  assumptions  like  finite 
cut-off  of  long  range  interactions.  With  the  exception  of  Xc-Pt(l  1 1)  system  (10)  there 
are  no  well  proven  or  tested  molecule-  or  ion-surface  potentials.  Consequently  we 
are  currently  unable  to  connect  directly  with  either  ex  situ  measurements 
e.g.,  I  on  Pt(lll)  or  Cl  on  Ag(100)  performed  by  transferring  emersed  electrodes 
into  U1IV  (11)  or  in  situ  surface  crystallographic  measurements  like  Sexafs  and  Xancs 
of  under  potentially  deposited  (upd)  metals  (12,  13). 

In  one  respect  the  model  used  here  was  simpler  than  the  set  of  assumptions  consti¬ 
tuting  the  standard  model  (2)  of  ionic  adsorption  on  metal  electrodes.  In  our  model 
all  species  were  attracted  to  the  electrodes  by  the  same  'wall'  potential.  This  potential 
and  the  static  electric  field  were  not  subject  to  the  0.82  run  cut-off  applied  to  all  the 
intermolecular  interactions. 
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The  S  i  2  parameters  describing  lithium  and  to  a  lesser  extent  fluoride  arc  suspect 
because  in  the  ease  of  lithium  ion  the  hydration  enthalpy  calculated  is  too  large 
compared  to  experiment  (14-16).  it  is  possible  that  the  adsorption  of  the  smallest  ions 
Lt  and  F  as  fully  hydrated  species  in  the  2x10  V/cm  Jicld  could  be  artifacts  due 
to  the  exaggerated  hydration  enthalpy.  However  this  argument  docs  not  affect  our 
general  conclusions  if  only  because  the  applied  static  field  used  is  at  the  upper  end 
of  fields  encountered  in  electrochemical  double  layers  and  if  more  realistic  model  pa¬ 
rameters  had  weaker  hydration  interactions  there  would  still  be  fields  where  the  small 
ions  would  adsorb  fully  hydrated.  In  weaker  fields  than  the  one  used  in  this  paper 
the  ions  relax  towards  film  center  (6).  We  will  discuss  the  effect  this  could  have  on 
the  interpretation  of  the  calculations  more  fully  in  a  later  section. 

Another  point  of  concern  is  the  electrical  discharge  of  contact  adsorbed  ions.  In  the 
older  papers  and  text  books  (2,  17)  halide  adsorption  is  treated  as  a  physisorption 
phenomenon.  However  there  is  exceptional  good  experimental  evidence  that  argues 
strongly  that  this  is  far  too  simplistic  a  picture.  After  contact  adsorption  a  large 
polarizable  ions  like  iodide  on  metals  Au  and  I’t  can  'discharge'  and  become 
chemisorbed  (11,  18-21).  Discharge  of  ions  and  the  formation  of  chemisorbed  species 
is  outside  the  scope  or  the  classical  molecular  dynamics  scheme  used  here, 
tains  its  full  charge. 

We  close  this  section  with  a  brief  review  of  recent  Ml)  calculations  in  the  general 
area  of  solid  -  aqueous  interfaces.  In  a  very  early  study,  the  Mainz  group  (22,  23), 
followed  the  evolution  of  eight  Li  I  and  200  water  molecules  between  Leonard  -  Jones 
walls  for  10  and  20  ps.  They  were  only  able  to  calculate  short  time  correlation 
functions  and  not  the  adsorption  process  which  required  a  longer  simulation  than  they 
were  able  to  do  at  the  time.  Rose  and  Benjamin  (24)  have  described  MI)  simulations 
of  a  single  Na  *  and  a  single  Cl  ion  in  a  slab  of  512  waters  (llcxiblc  SPC  model 


(25),  )  between  two  l’t(lOO)  surfaces.  The  geometry  and  potentials  were  those  devised 
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hy  llcinzingcr  and  coworkers  (26,  27).  l  or  l’t(IOU)  the  first  monolayer  of  water  was 
strongly  bound.  A  result  of  significance  was  the  calculation  of  the  Helmholtz,  free 
energy  of  adsorption  in  zero  electric  field.  The  anion  was  more  strongly  bound  than 
the  cation.  More  recently  Seitz- Bey wl  ct  al  (28)  have  studied  hydration  of  Li  and 
1  separately  in  three  layers  of  water  at  the  I’t  surface.  They  used  the  flexible  water 
model  of  Bopp  et  al  (29),  in  a  simulation  lasting  7.5  ps.  The  work  reported  here  is 
complementary  to  that  cited  and  described  above. 

There  is  also  continuing  interest  in  thin  films  of  water  and  electrolytes  in  connection 
with  a  possible  phase  transition  (30,  31),  vibrational  relaxation  dynamics,  hydrogen 
bond  formation  dynamics  (32,  33),  and  the  interaction  of  water  with  metals  treated 
as  jellium  half  spaces  (34). 

This  paper  is  arranged  as  follows.  Section  II  describes  the  models  used  for  ions  and 
water  and  details  of  computation.  Sections  Hi  through  VII  describe  the  results  for 
the  halide  ions,  and  the  last  section  summarizes  the  main  results  and  conclusions  for 
this  series  of  simulations. 


II.  THE  MODEL 
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In  the  calculations  described  here  we  use  the  water  model  originally  devised  bv 
Stillingcr  (7,  35),  and  commonly  referred  to  as  the  ST2  model.  Although  tins  model 
required  more  computer  time  because  it  contains  four  charged  and  one  uncharged 
point  particles,  this  draw  back  is  more  than  compensated  by  the  extensive  set  of 
Lcnnard-Joncs  parameters  for  use  with  all  the  alkali  metal  cations  and  halide  anions. 
These  parameters  were  developed  in  a  series  of  molecular  dynamics  calculations  bv 
llcinzingcr  and  Spohr  (26,  27). 

Interaction  Energy 

The  Coulomb  interaction  between  molecules  was  represented  as  sum  of  l/r  inter¬ 
actions  between  atomic  point  charges.  These  interactions  were  soften  for  small  mo¬ 
lecular  separation  by  a  switching  function  S(R)  .  This  switching  function  was 
introduced  by  Stillingcr  (7),  and  was  clearly  documented  by  Stcinhauscr  (35).  The 
short  range  part  of  the  intertnolccular  interaction  was  modeled  by  1. award-.} ones 
potential  between  the  atoms  of  each  molecule.  All  molecule-molecule  interactions 
(both  l/r  and  Lennard-.Ioncs  potentials)  were  cut-off  in  a  smooth  fashion  at  molecular 
separation  R  =  0.82  nm  by  a  truncation  function  7 (R)  .  The  atoms  of  each  molecule 
also  interacted  with  surfaces  at  z  ~  ±z0  .  Both  surfaces  were  treated  as  Hat  featureless 
plates  with  a  uniform  electric  charge  density  of  1-  7  and  —  n  on  the  —  z„  and  +  z„ 
plates  respectively.  This  gave  rise  to  a  uniform  electric  field,  J:  =  4nKo  in  the  positive 
z-dircction  where  the  electrostatic  coupling  constant  K  had  the  value  138.936 
KJ»nm/(moic«  e2)  in  the  units  used  in  this  calculation.  Non  Coulombic  interactions 
between  the  walls  and  all  the  atoms  were  represented  by  the  9-3  potential  (wall  po¬ 
tential)  introduced  by  Lee  ct  al  (36).  The  complete  interaction  energy  U  is, 
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where  /  and  _/  were  molecular  indices,  and,  a  and  /I  were  atomic  indices.  Here  A,  was 
the  set  of  II  atoms  of  molecule  /  ,  A„  was  the  distance  between  the  center  of  mass 
of  molecules  i  and  j,  and  The  symbol  r,p  was  the  distance  between  atoms  a  and  ft. 
For  small  A  the  Coulomb  energy  was  modified  by  the  switching  function 
S(R ,  A,,  A(/)  given  by, 
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The  values  of  At  and  Ay  were  dependent  on  the  types  of  the  molecular  species  that 
weic  interacting. 

As  mentioned  above  the  tails  of  Coulomb  and  l.cnnard-.loncs  pair  interactions  were 
cut  off  by  the  truncation  function  /(A).  The  form  of  /(A)  was  given  by, 


r  -  r; 

«?/-  Rr 


r<  it; 


R I  <  R  <  r'u 


R„<  R 


The  same  truncation  function  was  applied  to  all  molecular  interactions,  with 
R[  =  0.77V  Hi!  and  Rj,=  0.820  run.  The  integers  m  and  n  controlled  the  smoothness 
of  the  truncation  function  at  R]  and  R'u  respectively.  In  these  calculations  n  -  in  -  2 
which  insured  that  the  energy  had  smooth  first  spatial  derivatives. 

liquations  of  motion 

All  molecules,  in  this  study,  were  considered  as  a  rigid  collection  of  point  atoms.  That 
is  all  bond  lengths  and  angles  within  a  molecule  were  fixed.  1  o  evolve  a  collection 
of  these  molecules  a  quaternion  formulation  of  the  rigid  body  equations  of  motion 
was  used  (37).  The  center  of  mass  position  (R.)  and  velocity  (V,)  was  used  to  describe 
the  translational  degrees  of  freedoms  of  molecule  i.  The  orientational  motion  of  the 
molecule  was  described  by  the  quaternion  q,  =  (qf,  q],  q},  qf)  and  the  rotational  velocity 
(wf),  as  measured  in  the  body  frame  of  the  molecule.  The  one  exception  to  this  was 
for  monatomic  molecules,  in  which  ease  the  orientational  degrees  of  freedom  were  not 
needed. 

The  discussion  of  the  equations  of  motion  begin  by  considering  the  potential  energy 
IJ.  From  equation  1,  it  can  be  seen  that  the  potential  energy  can  be  treated  as  a 
scalar  function  of  tire  variables  (R,....,  R*;  r„),  where  N  is  the  number  of  molecules 

and  n  is  the  number  of  atoms  in  the  entire  system. 


U—  {7(Rt,...,  Rn)  rj,...,  r„) 


(d) 


Of  course  not  all  these  variables  are  independent.  However  (or  ihe  purposes  ol  the 


following  two  definitions  they  are  treater!  as  independent  variables 


f.  =  v;- 


f-vro. 


The  total  force  F,'  and  torque  r,,  can  he  expressed  in  terms  o|  f.  and  F,  as 


F"/  =  f,  +  Tt  —  ^<r„-  K,)xf„ 


The  translational  motion  of  tfic  molecule  is  described  by  the  first  order  ordinary  dif¬ 
ferential  equation. 


(/R. 

ur=v‘ 


<A;,  T 

“•IT  "  F<- 


For  the  rotational  motion  it  is  convenient  to  work  in  the  body  axis  of  the  molecule, 
where  the  moment  of  inertia  /,  is  diagonal  and  time  independent.  It  is  useful  to  define 

A 

the  operator  Q, 
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and  the  body  frame  rotational  force 


»-/?  b  b  r  wb  b\ 
F,  =  Tt  -  W,  X  (/jWj). 


Using  these  quantities  above,  the  following  first  order  ordinary  differential  equations 
can  he  written  to  describe  the  rotational  motion 


i. 


dt 


r.b 

"rfT  =  *'  ■ 


(IU) 


This  set  of  equations  conserve  the  total  energy  lor  time  independent  potentials  ((/)  . 
A  constant  temperature  ensemble  may  be  simulated  by  introducing  a  velocity  de¬ 
pendent  term  in  the  acceleration  terms  to  constrain  the  total  kinetic  energy.  The 
constant  temperature  equations  of  motion  are  written  as. 


dRt 

~7T 


v„ 


Mi 


d^L 

dt 


F/  -  yV; 


dt 


do/- 

dt 


i 


(ii) 


This  choice  of  y  ensures  that  the  total  kinetic  energy  of  the  systems  is  constant. 
Model  for  water  molecules  and  ions. 

In  the  ST2  model  the  water  molecule  consisted  of  a  central  oxygen  atom  (O  S  12) 
surrounded  by  two  hydrogen  atoms  (II  ST2)  and  two  massless  point  charges 
(PC_ST2)  in  a  rigid  tetrahedral  arrangement  (bond  angle  =  cos  '(!  iJT)  ).  The 
O  ST2/H  ST2  and  OST2/PCST2  bond  lengths  were  0.10  nm  and  0.08  nm  respec¬ 
tively.  The  oxygen  atom  was  the  only  Lennard-Joncs  atom’  in  ST2  water.  The 
hydrogen  and  point  charges  interacted  with  their  surroundings  (i.c.  atoms  and  sur¬ 
faces)  by  Coulomb  interactions  only.  The  ions  were  treated  as  non  polarizable 
Lennard-Jones  'atom'  with  point  mass  and  charge. 

Interaction  Parameters. 

The  atom-atom  and  atom-surface  interaction  parameters  arc  given  in  Table  1.  These 
parameters  arc  taken  from  work  of  llcin/.ingcr  and  Spohr  (26,  27).  The  combining 
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rules  for  unlike  species  were:  tAH  =  (>'AAr. aAB  -  Vi(aAA  +  anB).  The  switching  f  unc¬ 
tion  interval  ends  R'j  and  R‘{,  all  vanish  except  lor  S  I  2/S  1 2  pairs,  where 
/?£n-*n=0.20160  nm  and  0.3 1287  nm. 

Simulation  details 

The  systems  considered  had  either  216  or  432  molecules  and  ions.  In  the  main  set 
of  simulations  a  system  of  214  water  ST2  molecules,  one  Li  ion  and  one  halide  ion 
were  studied.  For  all  systems  the  molecules  were  confined  to  a  rectangular  simulation 
box  with  dimensions  1.862  nin  x  1.862  nm  x2.362  nm  {  x  x  y  x  z)  and  periodic- 
boundary  conditions  in  the  x  and  y  directions.  In  the  calculations  with  432  molecules 
the  simulation  box  had  the  same  x  and  y  dimensions  and  was  4.230  nm  thick  in  the 
z  direction.  Initially  the  molecules  were  randomly  disposed  on  a  cubic  lattice  with 
lattice  parameter  0.31  nm.  The  equations  of  motion  were  integrated  by  a  fourth  order 
multi-step  numerical  scheme  with  a  2  fs  time  step.  With  this  scheme  the  time  step 
to  time  step  rms  fluctuation  of  quaternions  length  and  cncrgy/tcmpcraturc  was  about 
0.0002%  and  0.0009%/0.002%  respectively.  Though  small,  these  changes  could  lead 
to  global  drift,  so  at  each  time  step  a  small  scaling  correction  was  made  to  the 
quaternions  and  velocities.  Also  the  global  center  of  mass  velocities  in  the  x  and  y 
directions  was  set  to  zero  at  each  time  step  by  shifting  the  molecular  translational 
velocities 

All  calculations  were  performed  on  an  IBM  RS/6000  model  540  operating  at  16 
Mflop.  For  the  simulations  involving  432  molecules  and  ions  a  100  ps  simulation 
with  2  fs  time  step  took  approximately  25  hours. 

Ion-water  pair  correlation  functions. 

The  ion  -  water  pair  correlation  functions  for  all  the  ions  were  calculated  in  the 
manner  described  in  the  first  paper  (1).  Figure  1  shows  the  result  for  all  the  ions 
studied  in  this  paper,  the  distribution  including  only  water  molecules  in  a  cone  with 
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half  angle  60°  pointing  into  the  bulk  with  axis  parallel  to  the  /  axis.  For  chloride, 
bromide,  and  iodide  the  second  peak  appeared  identical.  The  third  peak  in  the  figure 
contained  a  narrow  minimum  at  0.8  nm,  This  is  indicated  by  the  arrows  in  f  igure 
I.  This  feature  was  due  to  the  seam  in  the  water  distribution  due  to  the  truncation 
of  all  interactions  between  0.779  and  0.820  nm  (sec  Fqn.3). 

Cross  correlation  values. 

The  correlation  between  two  variable  v  and  tv  was  tested  by  evaluating  the  cross 
correlation  factor  /?, 


<  vw> 


<  v>  <  w> 


<  vv>  —  <  v>  2)(  <  ww>  —  <  w>  2) 


(12) 


where  <</>  >  denote  the  average  of  (f>  .  If  v  and  w  arc  completely  uncorrelatcd,  ft 
vanishes  and  if  v  and  w  are  highly  correlated,  the  magnitude  of  /?  will  be  close  to  one. 
The  cross  correlation  factor  was  used  to  explore  the  correlation  between  the  ion's 
position  and  hydration  number 

Residence  and  Exchange  Times 

It  was  of  interest  to  consider  the  mean  residence  time  of  a  water  molecule  in  a  ion's 
hydration  shell.  The  hydrations  shell's  inner  and  outer  radii  r ,  and  r0  respectively  are 
determined  by  the  minimums  of  the  ion-water  pair  correlations  function.  A  water 
molecule  is  consider  to  have  entered  the  hydration  shell  when  the  water-ion  separation 
r  satisfies  the  relationship,  r,(l  -F  S)  <  r  <  /•„(!  —  <5)  for  the  first  time.  The  parameter  ft 
is  a  small  number  that  filters  out  water  molecules  that  only  make  very  short  time 
excursions  into  the  shell.  For  this  calculation  <5  is  chosen  to  be  0.1.  A  water  mol¬ 
ecule  is  then  consider  to  have  left  the  hydration  shell  when  r<r, (!  —  <!>)  or 
r„(l  +  <  r  In  this  case  <5  filter  out  events  where  a  water  molecule  makes  a  short  time 

excursion  out  of  the  shell.  We  consider  the  residence  of  this  water  molecule  about 
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the  ion  to  be  the  time  between  these  two  events.  Then  mean  residence  time  is 
approximated  by  averaging  this  time  over  all  water  molecules  than  both  enter  and 
leave  the  hydration  shell  during  the  simulation. 

A  related  measure  of  the  rate  of  change  of  the  hydration  shell,  is  the  mean  time  of 
change  of  the  hydration  shell.  This  is  simply  the  average  time  between  events  where 
the  set  of  water  molecules  that  constitute  the  hydration  shell  changes.  In  this  paper 
a  hydration  shell  is  consider  to  have  changed  if  it  membership  differs  from  the  original 
for  more  than  Ips  A  stable  hydration  shell  is  considered  to  be  created  when  its 
membership  remains  unchanged  for  more  that  Ips.  1‘hc  mean  time  to  change  :.H  sub 
C  is  defined  as  the  time  to  change  of  a  stable  hydration  shell  averaged  over  all  stable 
shells  created  in  the  simulation. 
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III. RESULTS  FOF  LiF 

Figure  2  shows  the  time  evolution  of  the  z  coordinate  (see  left  side  vertical  scale)  of 
the  lithium  ion  and  fluoride  ion  for  a  time  period  lasting  1000  ps,  in  a  cell  also  con¬ 
taining  214  water  molecules.  The  broken  lines  at  |/|  =  0.9.14  nm  at  top  and  bottom 
of  the  figure  indicate  where  the  hard  wall  part  of  the  surface  potential  begins  i.c., 
where  V(z)  =  0.  The  configurations  used  in  constructing  this  figure  were  I  ps  apart. 
In  the  starting  configuration  both  ions  were  located  near  the  cathode  corresponding 
to  time  t  and  position  z  given  by  (t,  z)  =  (0,  -0.5).  In  f  igure  2  this  corresponds  to 
a  point  near  the  bottom  left  edge.  Under  the  influence  of  the  electric  field  the  fluoride 
ion  migrated  the  greatest  distance,  traversing  most  of  the  water  layer  to  reach  the 
anode  corresponding  to  the  broken  line  at  z  =  0.914  nm  at  the  top  of  the  figure. 
It  took  approximately  100  ps  for  the  fluoride  ion  to  traverse  the  film.  It  is  not  un¬ 
reasonable  to  assume  that  during  this  time  the  solvent  underwent  considerable  mixing, 
and  that  sufficient  time  lapsed  sufficient  time  for  the  whole  system  to  settle  into  the 
part  of  phase  space  describing  equilibrium  of  the  ions  and  associated  waters.  Notice 
that  in  the  applied  field  the  path  of  the  z  coordinate  of  the  fluoride  ion  for  the  first 
100  ps  was  not  a  monotonically  increasing  function  of  time.  The  initial  movement 
was  best  described  as  a  steady  drift  with  superimposed  high  frequency  jitter  most  likely 
due  to  interactions  of  the  ion's  primary  solvent  shell  with  surrounding  solvent.  This 
type  of  motion  has  been  frequently  seen  in  our  studies  of  hydrated  ions,  and  may 
well  be  typical  of  ions  pushing  through  a  ll-bondcd  network.  For  times  greater  than 
100  ps  the  time  dependence  of  the  z  coordinate  of  botli  ions  looks  qualitatively  the 
same.  The  range  of  positions  for  the  ion  was  about  0.4  nm,  not  much  more  than  the 
diameter  of  a  water  molecule  (a  =  0.11  nm).  flic  time  averaged  position  of  the 
fluoride  ion  was  -0.668  nm,  about  one  water  diameter  from  the  hard  wall  region  (|z| 
=  0.934  nm). 
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Also  shown  in  Figure  2  as  a  function  of  time  for  each  ion  are  the  numbers  of  water 
molecules  in  the  first  and  second  solvation  shells.  The  data  was  sampled  every  5  ps 
not  every  1  ps  as  for  the  z  coordinates.  The  radii  of  the  solvation  shells  were  chosen 
to  be  at  the  minima  in  the  calculated  ion-water  pair  correlation  functions.  1'or  lithium 
these  were  at  r  =  0.278.5  nm  (first  shell)  atid  r  =  0.500  nm  (second  shell).  Tor 
fluoride  these  were  at  r  =  0.3405  nm  (first  shell)  and  r  =  0.5527  nm  (second  shell). 

It  can  be  seen  that  for  lithium  the  first  shell  contained  six  waters  with  very  infrequent 
excursions  to  five  and  seven.  There  was  one  event  that  lasted  more  than  five 
picoseconds  (see  bottom  right  in  Figure  2).  The  time  averaged  first  hydration  shell 
number  was  n  =  6.01  .  This  is  consistent  with  the  calculation  of  Szasz  et  ai  (38), 
who  reported  the  same  value  for  bulk  electrolyte,  and  the  neutron  scattering  exper¬ 
iments  of  Newsome  et  al  (39),  and  I  Jewish  et  a!  (40).  The  mean  time  to  change  Tc 
of  the  first  solvation  shell  was  between  about  25  ps  and  35  ps.  This  was  in  the  range 
of  all  the  other  calculated  values  for  the  lithium  ion  with  different  anions.  Table  II 
summarizes  the  hydration  water  residence  times  for  all  the  ions.  In  the  case  of  lithium 
the  residence  time  was  about  100  ps  for  the  first  shell.  This  was  about  an  order  of 
magnitude  larger  than  for  the  second  shell,  and  approximately  consistent  with  the 
mean  time  between  changes  (25  x  6  vs.  106  ps).  These  numbers  arc  consistent  with 
those  reported  for  ions  in  a  bulk  environment  in  the  interesting  paper  of  Impey  et 
al  (15).  Roughly  the  residence  times  on  LiH  ions  in  different  LiX  solutions  were  the 
same,  with  the  larger  values  associated  with  greater  statistical  uncertainty.  It  was  also 
three  times  greater  than  residence  time  for  water  in  the  first  shell  of  fluoride  ion.  Tor 
the  first  shell  there  did  not  appear  to  be  much  correlation  between  position  and  hy¬ 
dration  number  for  the  first  shell.  The  calculated  correlation  coefficient  was  0.1.  In 
the  second  shell  the  hydration  number  fluctuated  between  seven  and  sixteen,  with  a 
calculated  average  value  of  11.4,  displayed  a  strong  correlation  with  distance  from  the 
electrode.  The  correlation  coefficient  of  0.5  supported  the  visual  impression.  The 


mean  time  to  change  of  the  second  shell  was  a  not  much  larger  than  the  sample  m 
terval  and  so  not  reliable,  except  to  say  that  it  was  of  the  order  of  1  ps.  Hie  picture 
that  emerged  was  that  the  lithium  ion  and  it's  first  shell  he  considered  as  a  rather  still 
or  even  as  a  rigid  object,  surrounded  by  a  compliant  second  shell. 

For  fluoride  the  first  shell  contained  seven  waters  with  frequent  excursions  to  six  and 
seven,  and  a  few  to  five.  The  average  first  shell  hydration  number  was  6.8.  The  mean 
time  between  changes  of  the  first  shell  was  3  ps.  This  was  an  order  of  magnitude 
smaller  than  for  lithium.  The  residence  time  (  see  Table  II)  for  hydration  water  was 
about  25  ps,  in  the  216  molecule  system,  and  a  little  longer  (33  ps)  in  the  larger 
system.  The  differences  arc  not  large  enough  to  be  considered  significant.  There  was 
no  apparent  correlation  between  first  shell  hydration  number  and  position,  and  this 
was  reflected  in  the  calculated  correlation  of  0.18.  In  the  second  shell  the  hydration 
number  fluctuated  between  ten  and  twenty  two,  with  an  average  value  of  15.4.  There 
was  a  correlation  with  distance  from  the  electrode  about  the  same  as  in  the  second 
shell  of  lithium.  The  correlation  coefficient  was  0.57. 

Figure  3  shows  the  density  profiles  for  lithium  and  fluoride  ions  (normalized  to  unity) 
and  water  (normalized  to  214)  across  the  gap.  Only  configurations  with  time  greater 
than  100  ps  were  used  in  accumulating  this  density  profile.  The  figure  also  shows 
schematically  a  cartoon  of  the  relative  positions  of  the  ions  and  water  molecules  in 
the  first  hydration  shell  closest  to  the  surface,  (dearly  in  the  density  peak'  position 
both  ions  are  hydrated.  Not  surprisingly  the  cation  has  a  peak  closer  to  the  electrode 
since  it  is  the  ion  with  the  smaller  radius.  The  lithium  peak  position  was  at  -0.709 
nm  and  the  fluoride  was  at  0.661  nm.  The  average  positions  were  -0.668  nin  (Li4 ) 
and  0.628  nm  (F).  The  difference  0.05  nm  should  be  compared  to  the  difference  in 
Lennard-Jones  radii  which  is  0.08  nm.  The  range  of  positions  of  both  ions  is  very 
similar,  a  little  less  than  0.5  nm  in  each  ease.  The  hydrated  ions  have  the  density 
profile  of  a  particle  trapped  in  a  well  undergoing  Brownian  motion.  For  convenience 
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the  main  features  of  the  densities  for  ions  and  water  arc  collected  into  Tables  111  and 
IV. 

The  water  density  profile  lias  interesting  features.  The  larger  peak  at  -0.803  nm  was 
due  to  water  in  the  first  layer  of  the  sample.  When  structure  in  the  water  distribution 
is  compared  with  the  ions  it  is  clear  that  the  larger  size  of  the  anion  is  responsible 
for  the  shallow  minimum  located  near  0.66  nm.  This  conclusion  was  further  sub¬ 
stantiated  by  the  increase  in  size  of  this  feature  with  size  of  the  anion.  This  wiil  be 
seen  clearly  in  later  figures  of  density  profiles.  Calculations  of  the  density  profiles  for 
water  center  of  mass  and  the  hydrogen  atoms  were  also  performed.  The  most  inter¬ 
esting  feature  of  these  calculations  was  a  small  but  distinct  peak  in  the  hydrogen 
distribution  near  the  cathode  at  -0.992  nm.  This  was  due  to  protons  oriented  towards 
the  electrode. 

Figure  4  shows  density  profiles  calculated  with  the  wider  box  (4.230  nm)  containing 
twice  as  many  molecules  and  ions.  The  hard  wall  was  located  at  z.  =  ±  1.868  nm. 
Notice  that  the  density  profiles  are  narrower  (smaller  full  width  at  half  maximum). 
The  reason  for  this  was  traced  to  smaller  fluctuations  in  the  positions  of  the  ions 
because  of  weaker  transcell  interactions.  The  plots  of  z  vs.  t  (time)  did  not  show  the 
same  range  of  movement  seen  in  Figure  2.  Also  the  ions  were  more  strongly  adsorbed 
as  judged  by  their  peak  positions,  which  were  each  closer  to  their  respective  electrodes. 
The  peak  positions  relative  to  the  hard  wall  are  listed  in  Table  III.  For  fluoride  the 
distance  from  peak  to  hard  wall  was  0.13ntn  compared  to  0.27  nm  for  the  thinner  film 
with  216  molecules  and  ions.  The  cartoon  for  the  fluoride  ion  shows  in  to  be  closer 
to  the  electrode  in  Figure  4  compared  to  Figure  3.  There  was  also  a  small  shift  to 
the  electrode  Tor  the  cation,  but  it  was  smaller  due  to  stronger  binding  of  the  waters 
of  hydration.  With  a  cut  off  radius  of  0.82  nm  the  ions  did  not  couple  directly  by 
Coulomb  interactions.  There  was  no  overlap  of  the  range  of  the  interactions  that 
would  be  possible  at  times  in  the  2.362  nm  thick  sample. 
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Therc  is  some  concern  that  the  transccll  ion  attraction  in  the  thin  films  may  he  an 
artifact  due  to  the  truncation  of  the  Coulomb  potentials  of  the  sort  recently  described 
by  Hadcr  and  Chandler.  (41).  Bader  and  Chandler  identified  an  artificial,  truncation 
induced  attraction  between  ion  of  like  charge,  however  argue  that  for  ions  of  unlike 
charge  such  artifacts  should  be  reduced.  This  and  the  fact  that  the  transccll  attraction 
is  decreasing  with  film  thickness  argue  that  we  arc  not  seeing  such  truncation  induce 
effects. 

Figure  4  also  showed  that  the  water  at  the  interface  was  quieter,  as  measured  by  the 
peaks  and  troughs  in  the  water  density  profile.  Major  feature  positions  arc  recorded 
in  Table  IV.  Note  first  that  the  density  in  the  middle  of  the  films  was  about  the  same. 
The  peaks  corresponding  to  surface  water  were  higher  at  both  surfaces  relative  to  bulk 
for  the  432  molecule  system  compared  to  the  216  system.  Corroborating  this  notion 
was  the  deeper  water  exclusion  trough  near  the  fluoride  ion,  implying  a  higher  degree 
of  time  averaged  organization  around  the  anion. 
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V.  RESULTS  FOR  LiCI 

Figure  5  shows  the  time  evolution  of  the  z  coordinate  for  the  lithium  ion  and  chloride 
ion  in  a  216  molecule  and  ion  simulation.  Shown  also  is  the  time  dependence  of  the 
number  of  water  molecules  in  the  first  and  second  solvation  shells  for  chloride. 
Lithium  is  not  shown  because  it  was  essentially  identical  to  the  result  already  plotted 
in  Figure  2.  The  starting  configuration  was  the  same  as  previously  described  for  LiF. 
Note  first  that  the  time  dependence  of  the  lithium  ion  was  qualitatively  the  same  as 
in  the  LiF  simulation  and  is  not  discussed  further  here.  The  time  dependence  of  the 
z  coordinate  for  chloride  was  different  to  fluoride.  During  the  first  100  ps  of  the  LiCI 
simulation  hydrated  Cl  moved  faster  under  the  influence  of  the  electric  field  and  ar¬ 
rived  at  the  anode  sooner  than  fluoride  did  under  similar  conditions.  This  observation 
is  significant  since  it  is  in  qualitative  agreement  with  experimentally  measured  trans¬ 
port  numbers  in  dilute  solutions  (42-44).  This  difference  was  due  to  the  larger  effec¬ 
tive  radius  of  hydrated  fluoride  which  binds  its  solvation  shells  more  tightly  than  the 
chloride  ion.  In  passing  we  note  that  the  chloride  ion  z  coordinate  also  showed  jitter 
similar  to  that  already  described  for  fluoride  in  Figure  2. 

For  times  greater  than  100  ps  the  time  dependence  of  the  z  coordinate  of  chloride 
looks  qualitatively  different  from  the  strongly  hydrated  ions  fluoride  and  lithium.  First 
note  that  the  chloride  nucleus  occasionally  touches  and  penetrates  the  hard  wall, 
which  in  zero  electric  field  occurs  at  |zj  =  0.934  nm.  This  means  that  chloride  spends 
some  time  contact  adsorbed,  in  contrast  to  fluoride  which  never  made  direct  physical 
contact  with  the  electrode. 

Inside  the  box  in  Figure  5  the  data  for  the  first  and  second  solvation  shells  of  chloride 
was  sampled  every  5  ps  not  every  1  ps  as  for  the  z  coordinates.  The  radii  of  the 
solvation  shells  were  chosen  to  be  at  the  minima  in  the  calculated  chloride  ion-water 
pair  correlation  functions,  f  or  chloride  these  were  at  r  -  0.3758  nm  (first  shell)  and 
r  =  0.6558  nm  (second  shell). 
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For  chloride  the  first  shell  lluctuatcd  mostly  in  the  range  live  to  seven  with  excursions 
to  three  and  occasionally  as  high  as  ten  (not  shown  in  the  time  spanned  in  figure  5). 
The  average  hydration  numbers  for  the  first  and  second  hydration  shells  averaged  over 
900  ps  were  6.1  and  20.8  respectively.  The  first  number  is  consistent  with  the  pio¬ 
neering  studies  of  Dietz  et  al  (45),  who  calculated  6.5  ±0.5  for  the  chloride  ion  in 
bulk  electrolyte.  There  was  an  apparent  correlation  between  first  shell  hydration 
number  and  position  when  the  ion  was  close  to  the  electrode.  The  calculated  corre¬ 
lation  was  0.51.  For  times  greater  than  100  ps  the  second  shell  hydration  number 
fluctuated  between  twelve  and  thirty  seven.  The  averaged  hydration  number  was  20. 8 
for  the  second  shell.  Again  there  was  apparent  correlation  with  distance  from  the 
electrode  and  the  calculated  correlation  was  0.81,  higher  than  for  the  shells  of  smaller 
ions.  The  residence  time  for  water  on  chloride  was  about  15  ps  for  the  large  system 
and  a  little  smaller  for  the  small  system  (see  Table  II).  It  was  not  dependent  on  shell 
indicating  that  the  shells  are  mixing  in  comparable  fashion  consistent  with  water  being 
weakly  bound.  This  was  the  common  description  for  all  except  the  two  smallest  ions. 

Figure  6  shows  the  z  dependent  density  profiles  for  lithium  and  chloride  ions  (nor¬ 
malized  to  unity)  and  water  (normalized  to  214)  across  the  film.  Only  configurations 
with  times  greater  than  100  ps  were  used  in  accumulating  these  density  profiles.  The 
total  simulation  time  was  actually  1100  ps  more  than  shown  in  the  figure,  so  that  for 
chloride  statistics  were  collected  for  a  total  duration  of  1000  ps.  The  main  peak 
(average)  of  lithium  ion  occurred  at  -0.732  nm  (-0.679  nm)  with  closest  approach  at 
-0.827  nm  and  furthest  position  at  z  =  -0.283  nm.  Values  arc  tabulated  in  fable  III. 
The  lithium  distribution  resembles  that  for  LiF  quite  closely. 

For  chloride  the  main  peak  (average)  was  at  0.709  nm  (0.747  nm),  closest  at  0.968 
nm  and  farthest  point  at  z  =  0.449  nm.  The  range  of  displacements  of  the  chloride 
ion  was  about  0.5  nm.  This  was  greater  than  shown  in  Figure  4  which  only  showed 
the  first  nanosecond,  a  large  excursion  away  from  the  electrode  occurred  between  1000 
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ps  and  1100  ps.  (dearly  the  statistics  were  not  sufficient  to  give  a  smooth  distribution, 
and  it  is  possible  that  the  density  is  bimodal  t or  a  gap  of  2.362  nm.  If  the  density 
profile  is  bimodal  then  a  calculation  of  the  potential  of  mean  force  (46)  will  show  two 
minima.  This  latter  conclusion  was  also  suggested  by  the  rapidity  of  the  transition 
from  contact  to  non-contact  adsorbed  configuration  (see  Figure  4).  It  was  considered 
not  worthwhile  to  pursue  this  aspect  further  with  the  model  used  in  this  paper.  The 
figure  also  shows  schematically  a  cartoon  of  the  relative  positions  of  the  ions  and 
water  molecules  in  the  first  hydration  shell  closest  to  the  surface.  Only  when  the 
chloride  ion  is  furthest  from  the  electrode  at  z  =  0.5  nm,  is  it  possible  for  the  ion 
to  be  surrounded  by  a  complete  first  solvation  shell.  1'hc  position  of  closest  approach 
is  slightly  greater  than  0.934  nm  because  the  attractive  electric  field  acts  to  pull  the 
ion  into  the  repulsive  region.  Also  shown  in  Figure  6  is  the  water  density  distribution. 
Chloride  being  a  larger  anion  than  fluoride  gives  rise  to  a  deeper  minimum  located 
near  0.66  nm  the  position  of  the  main  peak  (at  7.  =  0.71  nm)  in  the  ion  density. 
Other  features  in  common  with  the  fluoride  simulations  were  the  sharpness  of  the 
surface  water  peaks  (  at  |z(  =  0.898  nm  )  and  a  separate  peak  in  the  calculated  proton 
distribution  at  z  =  -0.992  nm  near  the  anode  (not  shown  in  Figure  6). 

Figure  7  displays  the  density  profile  results  for  the  thicker  4.230  nm  film.  Note  the 
dramatic  change  in  distributions  for  both  ions.  Hach  was  found  to  be  sharper  and 
occupying  a  narrower  range  of  Az.  The  full  widths  at  half  maximum  were  comparable 
to  those  found  for  all  the  other  thicker  systems.  The  cartoon  of  hydration  around 
the  chloride  ligand  in  Figure  7  showed  that  the  anion  was  more  strongly  adsorbed 
than  in  the  2.3  nm  film.  The  half  widths  and  position  shifts  arc  consistent  with  the 
view  that  transccll  coupling  mediated  by  water  polarization  is  responsible  for  the  much 
broader  distributions  seen  in  the  2.36  nm  films.  The  chloride  distribution  remained 
wider  than  any  other  ion  in  similar  cells,  suggestive  of  the  behavior  seen  in  the  thinner 
film.  Chloride  behaved  intermediate  between  adsorbed  as  a  l  ydrated  ion  and  adsorbed 
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in  contact.  The  water  exclusion  zone  near  /  -  1.6  nm  resembled  that  seen  already 
for  the  fluoride  ion  in  Figure  4.  All  the  water  density  profiles  were  remarkably  similar, 
a  fact  reflected  in  part  by  the  similarity  in  the  values  in  Table  IV. 


VI  RFSl 1  IS  FOR  I  iiir 
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Figure  8  displays  the  time  dependence  ol  the  /  coordinate  ol  the  lithium  ion  and  a 
bromide  ion  during  a  simulation  lasting  M>0()  ps.  Unlike  all  the  simulations  presented 
so  far  this  one  started  with  an  equilibrated  system  in  which  the  ions  were  adsorbed 
on  the  electrodes.  It  was  (bund  that  starting  with  the  cubic  lattice  lead  to  the  for¬ 
mation  of  an  ion  pair.  The  procedure  used  for  bromide  was  to  start  with  the  6  x  6 
x  6  cubic  array  as  before  but  in  an  electric  field  ten  times  stronger,  flic  simulation 
was  run  in  high  field  for  10  or  20  ps  until  both  ions  were  in  contact  with  the 
oppositely  charged  wall.  In  fields  of  1.5  x  10  V/ cm  and  greater  even  the  cation  was 
contact  adsorbed.  With  the  ions  separated  the  field  was  then  reset  at  2xl07  V/cm. 
Dashed  horizontal  lines  define  the  hard  wall  region  of  the  wall  potential  starting  at 
|z|  =  0.934  nm.  This  is  effectively  the  boundary  surface  of  the  electrode.  Note  that 
the  bromide  ion  frequently  touched  the  boundary,  and  so  in  this  model  is  a  moderate 
contact  adsorber,  more  strongly  bound  than  chloride  in  Figure  4,  but  not  as  strongly 
bound  as  iodide  (vide  ultra). 

The  time  dependence  of  the  z  coordinate  of  the  lithium  ion  was  a  little  different  to 
that  found  for  FiCl  and  FiBr.  The  range  of  displacements  is  clearly  smaller,  which 
results  in  the  distribution  being  narrower.  We  return  to  this  subject  later. 

Shown  in  Figure  8  as  a  function  of  time  for  each  ion  arc  the  numbers  of  water 
molecules  in  the  first  and  second  solvation  shells  of  the  bromide  ion.  The  data  was 
sampled  every  5  ps  not  every  1  ps  as  for  the  z  coordinates.  The  radii  of  the  solvation 
shells  were  chosen  to  be  at  the  minima  i.i  the  calculated  ion-water  pair  correlation 
functions.  For  bromide  these  were  at  r  =  0,4024  nm  (first  shell)  and  r  =  0.6588  nm 
(second  shell),  same  as  for  chloride. 

l  or  bromide  the  contents  of  the  first  shell  fluctuated  between  live  and  six  waters  wi'h 
infrequent  excursions  to  four  and  occasional  jumps  to  values  as  high  as  nine.  The 
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averagc  value  was  5.6  lor  the  lirst  shell.  The  higher  values  were  only  weakly  corre¬ 
lated  with  positions  further  from  the  electrode.  The  correlation  coellicicnt  lor  the  first 
shell  was  calculated  to  be  0.37.  In  the  second  shell  the  hydration  number  fluctuated 
between  ten  and  twenty  six.  The  hydration  average  value  was  16.8  and  the  correlation 
coefficient  was  0.86,  implying  a  strong  correlation  with  distance  from  the  electrode. 
Table  II  shows  the  residence  times  follow  the  pattern  established  for  chloride. 

figure  9  shows  the  distribution  of  Li',  Br ,  and  214  water  molecules  in  a  strong 
electric  field  (  2xl0?  V/ cm)  across  the  gap  between  the  charged  plates.  The  wall  po¬ 
tential  is  shown  for  reference.  Circles  have  the  l.cnnard-.! ones  radii  for  the  ions  and 
water.  The  configurations  were  collected  in  the  time  interval  100  ps  <  time  <  1000 
ps.  The  range  of  bromide  positions  was:  farthest  z  =  0.661  nm,  closest  z  =  0.968 
nm,  and  peak  (average)  position  at  approximately  0.874  nm  (-0.861  nm).  At  the 
closest  approach  the  negative  ion  lay  inside  the  repulsive  regime  of  the  wall  potential. 
In  this  simulation  the  anion  engages  in  moderate  contact  adsorption,  stronger  than 
chloride  but  not  as  strong  as  in  the  ease  of  iodide.  Note  the  again  the  presence  of 
the  shallow  minimum  in  the  water  distribution  near  the  average  position  of  the  anion. 
This  minimum  is  slightly  more  pronounced  than  in  the  ease  of  chloride.  It  is  due  to 
exclusion  of  water  molecules  from  the  interface  by  the  larger  volume  of  the  ion. 

The  results  of  simulations  with  a  twelve  (4.230  nm)  water  layer  system  arc  shown  in 
Figure  10.  The  trend  to  contact  adsorption  with  larger  ionic  size  was  continued. 
Again  there  was  a  narrow  lithium  distribution  characteristic  of  the  thicker  film.  What 
was  of  interest  in  comparing  the  density  profiles  in  Figures  9  and  10,  was  their  simi¬ 
larity.  In  all  the  other  LiX  systems  the  2.36  and  4.23  nm  films  show  how  the  transccll 
interaefion  drops  in  going  to  the  thicker  cell. 

The  transccll  interactions  between  Li '  and  Br  appeared  weaker  in  thin  cell  compared 
say  to  LiF  and  LiCI.  As  with  Li  I ,  to  be  discussed  in  the  next  section,  it  appeared 
that  as  contact  adsorption  increased  then  the  more  the  anion  decoupled  from  the  bulk 
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watcr  and  the  smaller  the  perturbation  on  the  cation  due  to  movement  of  the  anion 
perpendicular  to  the  surface.  As  seen  from  the  Figure  10  and  the  positions  tabulated 
in  Table  III  the  lithium  cation  peak  at  z  =  -1.692  nm  had  moved  a  little  closer  to 
the  hard  wall,  and  the  density  profile  was  somewhat  sharper.  Indeed  the  lithium 
position  was  close  to  that  for  a  lithium  ion  without  an  anion  but  trapped  against  the 
cathode  by  the  static  field  used  in  all  these  simulations  (6). 
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VII.  RESULTS  FOR  Lil 

The  ease  of  lithium  iodide  solutions  at  constant  energy  were  described  in  the  first  (1) 
paper.  New  results  are  presented  here  for  simulations  done  at  constant  temperature. 
Iodide  represents  the  extreme  case  of  contact  adsorption.  Experimentally  iodide  is 
known  to  contact  adsorb  on  metals  at  potentials  corresponding  to  thermodynamically 
stable  electric  double  layers,  in  contrast  to  lithium  ions  which  do  not.  lime  de¬ 
pendence  of  the  z  coordinate  of  the  lithium  ion  and  a  iodide  ion  during  a  simulation 
lasting  400  ps  is  shown  in  Figure  11.  The  sampling  interval  in  these  calculations  was 
0.5  ps,  smaller  than  in  all  the  previous  calculations  described  so  far  which  were  all 
at  1.0  ps.  The  starting  configuration  was  one  previously  obtained  in  the  earlier  study 
(1).  As  before  the  dashed  horizontal  lines  define  the  repulsive  region  of  the  wall  po¬ 
tential  starting  at  |z|  =  0.934  nm.  'This  is  effectively  the  boundary  surface  of  the 
electrode.  Note  that  the  iodide  ion  spends  almost  all  its  time  at  the  boundary,  and 
in  this  model  is  classified  as  a  strong  contact  adsorber.  It  makes  only  short  time 
excursions  away  from  the  surface,  and  is  rarely  found  more  than  0. 1  nm  removed  from 
contact. 

Also  shown  in  Figure  11  as  a  function  of  time  for  the  iodide  ion  are  the  numbers 
of  water  molecules  in  the  first  and  second  solvation  shells.  The  data  was  sampled 
every  2.5  ps  not  every  5  ps  as  for  the  hydration  number  plots  in  the  figures  for  LiF, 
LiCl,  and  LiBr.  The  radii  of  the  solvation  shells  taken  from  the  calculated  ion-water 
pair  correlation  functions  were  at  r  =  0.4289  nm  (first  shell)  and  r  =  0.6765  nm 
(second  shell). 

For  iodide  the  first  shell  fluctuated  between  two  and  eight,  with  four,  five  and  six 
being  the  most  frequent  occupation  number.  The  average  value  was  5.4  for  the  first 
shell  with  a  correlation  of  0.29  (weak,  corresponding  to  no  apparent  correlation  visible 
in  Figure  11)  between  first  shell  hydration  number  and  position.  Because  the  ion  was 
so  strongly  adsorbed  the  correlation  was  small.  In  the  second  shell  the  hydration 
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number  fluctuated  between  eleven  and  twenty  live.  The  average  value  was  16.8.  A 
weak  correlation  with  distance  from  the  electrode  was  apparent,  the  calculated  corre¬ 
lation  was  0.63. 

As  shown  in  Figure  12  the  essential  features  of  the  iodide  distribution  emerge  even 
after  a  simulation  running  for  400  ps.  The  distribution  of  the  Li'  ion  and  the  214 
water  molecules  looks  very  similar  to  LiF  and  MCI  solutions.  In  this  simulation  the 
anion  exhibits  strong  contact  adsorption.  The  anion  density  profile  was  sharply 
peaked  (average)  at  z  =  0.921  nm  (0.886  nm)  very  close  to  the  hard  wall.  The  point 
of  closest  approach  was  0.992  nm,  and  that  farthest  from  the  electrode  was  at  0.685 
nm.  The  water  distribution  in  Figure  9  has  more  noise  than  the  one  shown  in  Figure 
6  for  LiCl  solution  because  of  poorer  statistics  resulting  from  a  simulation  time  of 
400  ps  duration.  In  spite  of  this  the  shallow  minimum  in  the  water  distribution  due 
to  the  ion  displacing  water  was  clearly  visible  near  0.70  nm.  Table  IV  gives  the  salient 
features  of  all  the  water  distributions. 

In  Figure  13  the  results  of  the  thick  film  density  profile  arc  displayed.  As  before  the 
system  was  allowed  to  equilibrate  for  100  ps  before  averaging  over  configurations 
collected  during  the  next  900  ps.  The  iodide  density  distribution  was  the  sharpest, 
of  all  the  halide  profiles  for  all  the  systems  studied  here.  The  peak  value  of  nearly 
14  nm  1  in  the  thick  film  compared  to  1 1  nm  1  in  the  six  layer  film.  A  comparison 
of  water  density  profiles  showed  them  all  to  be  very  similar,  with  all  the  distribution 
differences  occurring  near  the  walls.  The  differences  that  did  exist  seemed  to  be  re¬ 
lated  to  the  degree  and  strength  of  hydration  of  the  adsorbed  ions.  Consistent  with 
weaker  coupling  across  this  thicker  film,  the  water  distribution  near  the  cathode  where 
the  Li  was  adsorbed  were  superimposabic  for  different  anions  in  systems  with  432 


molecules  and  ions. 


VIII.  RESULTS  FOR  WEAKER  FIELDS 


This  section  resumes  the  discussion  of  the  quality  of  water  and  ion  parameters  for 
surface  studies  begun  in  the  introduction. 

It  has  been  pointed  out  (8)  that  small  ion  parameters  over  estimate  the  of  hydration 
by  as  much  as  25%.  Consequently  in  the  field  wc  arc  using  2xl0?  V/cm  in  our  sim¬ 
ulations,  the  ability  of  small  ions  like  fluoride  F  and  lithium  Li  to  adsorb  with  intact 
primary  solvation  shell  is  exaggerated  compared  to  the  larger  ions.  The  magnitude 
of  this  field  is  at  the  upper  end  for  stable  electrochemical  double  layers  so  it  is  of 
interest  to  know  what  happens  in  weaker  fields. 

However  it  is  an  experimental  fact  that  small  ions  like  Li'  and  fluoride  1'  do  not 
chemisorb  on  to  metals  in  the  double  region  of  electrode  potentials,  and  this  phe¬ 
nomenon  is  attributed  to  their  strongly  held  primary  solvation  shells.  There  arc  even 
anions  (PF  BF  '  and  CIO  ")  that  adsorb  more  weakly  than  fluoride  (5),  in  that  they 
distrurb  the  interface  even  less  (lower  capacitance).  Consequently  although  the  ST2 
parameters  for  Li+  do  not  qualitatively  model  the  adsorption  process  at  the  value 
of  the  static  field  used  here,  the  scries  of  Ncinzingcr's  anion  parameters  does  model 
adsoption  changes  along  the  series.  We  do  not  expect  this  modelling  to  be  quanti¬ 
tative,  for  example  by  predicting  the  value  of  the  threshold  Held  above  which  the 
lithium  ion  to  be  contact  adsorbed. 

To  explore  the  static  field  dependence  wc  calculated  the  density  profile  for  a  single 
Li+  ion  with  215  water  molecules  and  a  single  1  ion  in  215  ST2  waters,  for  different 
values  of  the  applied  field.  All  other  conditions  in  these  simulations  were  the  same 
as  before. 

For  fields  about  seven  times  higher  (ca.  1.5x10  V/cm),  like  the  ones  used  in  the  cal¬ 
culations  of  Brodsky  and  Reinhardt  (30,  31),  the  lithium  ion  was  contact  adsorbed 
and  it  was  observed  that  all  the  waters  of  the  film  oriented  with  dipoles  ||/..  The 
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bchaviour  in  the  regime  of  weaker  fields  is  more  relevant  to  the  present  study.  As 
the  field  was  weakened  the  density  profile  of  Li  ion  shifted  towards  the  middle  of 

the  water  film,  in  zero  field  the  density  profile  p  (/.)  was  a  broad  bell-like  distrib- 

ion 

ution  with  the  tails  reaching  out  to  about  |z|  =  0.5nm.  This  result  is  consistent  with 
that  reported  by  Wilson  ct  al  (47)  for  a  single  sodium  ion  in  a  2.0  nm  thick  film  of 
water  with  two  liquid  -  vapour  boundaries,  and  can  be  interpreted  as  the  ions  mi¬ 
grating  to  maximize  their  polarization  energy  in  the  solvent.  Relaxation  times  were 
found  to  be  longer  than  1  ns,  because  of  the  large  configuration  space  explored. 
Consequently  it  was  not  possible  to  get  good  statistics  by  accumulating  configurations 

with  brute  force  dynamics.  This  problem  could  be  better  tackled  using  umbrella 

% 

sampling  (24)  or  some  scheme  to  simplify  the  long  time  dynamics  like  those  used  to 
studied  polypeptides  (48). 

In  the  case  of  iodide  the  peak  in  the  density  profile  p  (z)  for  a  given  field  always 

ion 

occurred  at  a  point  closer  to  the  charged  wall  than  was  the  case  for  lithium  ion.  This 
is  consistent  with  the  previous  calculations  for  I  ion,  and  suggests  that  the  estab¬ 
lishment  of  the  solvent  structure  around  this  large  anion  occurs  gradually  as  the  field 
is  weakened  and  the  ion  pulls  pulls  back  from  the  wall. 

What  can  be  concluded  from  these  results?  First  because  of  differences  in  ionic  radii 
the  smaller  ions  will  be  more  strongly  hydrated  than  larger  ions.  As  the  strength  of 
the  field  is  increased  ions  with  more  strongly  bound  solvation  shells  will  tend  to  spend 
more  time  contact  adsorbed.  Conversely,  in  the  absence  of  a  strong  universally 
binding  interaction  such  as  an  image  interaction,  the  adsorbed  ions  will  shift  away 
from  the  electrode  as  the  field  is  weakened.  An  exception  will  occur  for  large 
hydrophobic  ions  which  might  remain  adsorbed  in  a  configuration  that  maximized 
solvent  entropy  (it  is  well  established  that  large  neutral  organics  eg.,  benzene,  arc 
adsorbed  most  strongly  at  the  potential  of  zero  charge)  (49). 
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IX.  SUMMARY 

The  simulations  described  here  have  revealed  some  of  the  static  and  dynamical  com¬ 
plexities  of  the  adsorption  of  hydrated  ions  onto  charged  surfaces.  The  main  features 
were  as  follows,  i)  The  small  ions  Li  *  and  F  did  not  make  contact  with  electrode 
in  fields  of  2xlO;V/cm  because  of  their  strongly  bound  first  solvation  shells.  There 
was  a  balance  between  electric  and  Lcnnard-Joncs  attraction  and  thermal  fluctuations. 
In  the  thin  film  polarization  coupling  mediated  by  water  molecules  increased  the  size 
of  the  fluctuations  experienced  by  the  ions,  ii)  In  the  model  used  here  chloride  was 
the  first  ion  to  exhibit  contact  adsorption.  In  the  2  nm  film  it  appeared  to  have  a 
bimodal  distribution,  iii)  The  presence  of  adsorbed  ions  implies  the  presence  of  lo¬ 
calized  water  held  in  place  by  the  field  of  the  ion  and  it's  associated  hydrogen  bond 
network.  The  presence  of  large  interfacial  ions  also  means  there  is  a  depletion  zone 
where  fewer  waters  are  located,  iv)  The  time  dependence  of  the  first  shell  hydration 
numbers  for  lithium  and  fluoride  corroborate  the  picture.  The  residence  time  of  a 
water  molecule  on  lithium  ion  was  a  factor  of  four  longer  compared  to  fluoride.  The 
first  solvation  shell  of  chloride  and  the  bigger  ions  were  found  to  be  all  very  dynamic 
with  residence  times  around  10  ps  according  to  the  model  used  here,  exchange  of 
water  between  the  first  shell  and  the  bulk  was  less  than  about  1  ps. 

Finally  we  point  out  that  their  arc  significant  physical  effects  not  modeled  in  the 
present  study.  In  particular  rigid  non  polarizable  models  exhibit  dielectric  saturation. 
To  more  accurately  account  for  relaxation  around  small  ion  core  one  clearly  needs 
to  enhance  the  models  in  significant  ways  outside  the  scope  of  the  present  study.  In 
the  extremely  high  fields  near  small  ion  cores  water  molecules  will  distort  electronically 
and  geometrically.  To  account  for  these  effects  commonly  adopted  enhancements 
include  replacement  of  the  repulsive  r  12  core  potentials  with  softer  ones  fashioned 
after  exp(-ar),  or  better  still  a  sum  of  exponentials  obtained  from  a  quantum  chemistry 
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calculation,  polarizable  ions  and  neutral  water  molecules,  and  flexible  water  molecules. 
Work  on  incorporating  these  effects  in  useable  models  is  in  progress  in  our  laboratory. 
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'I'ablc  I  The  interaction  parameters  (q,r.,a.  A,  H)  and  mass  (»i)  (or  all  atoms  used 
in  the  simulations.  where  </„  =  <’  .  r.„  —  1  KJjmoie  ,  n„  I  mu  , 
A„  —  17.447  x  10  *KJ  •  nnflmolc  ,  B„  —  76.144  x  iO  iKJ  •  nni'jmole  .  and  m0  —  l  A  MU. 
The  e  and  a  for  unlike  atom  pairs  is  formed  from  the  combination  rules, 
£ab  —  \[(zaaZrb)  and  =  (csaa  + 


<tlq° 

ajo„ 

rnjm0 

AIA„ 

n\n 

OST2 

0.0000 

0.316 

0.310 

16.0 

1 

i 

1 1  ST2 

0.2357 

0.000 

0.000 

1 .000 

0 

0 

PCST2 

-0.2357 

0.000 

0.0 

0.0 

0 

0 

Li 

1 .0000 

0.149 

0.237 

6.9 

1 

1 

F 

-  1.0000 

0.050 

0.400 

19.0 

1 

f 

t 

Cl 

-  1.0000 

0.168 

0.486 

35.5 

1 

1 

Br 

-  1.0000 

0.270 

0.504 

79.9 

1 

1 

I 

-  1.0000 

0.408 

0.540 

129.9 

1 

1 

Tabic  11.  Summary  of  hydration  water  residence  times  wit  it  the  otic  sigma  errors  lot 
I.iX  (X  =  F,CI,Br,I)  in  water  films  with  216  molecules  (2.362nm)  and  432  molecules 
(4.320nm). 


LiF 

LiC  1 

LiBr 

Li  I 

2.362  nm 

film 

with  216 

molecules 

Li+  ion. 

1st 

shell 

122+25 

168+43 

115+43 

Li+  ion, 

2nd 

shell 

12+0.5 

12+0.4 

12+0.5 

X"  ion. 

1st 

shell 

25+1.8 

12+0.6 

11+0.6 

X  ion , 

2nd 

shell 

13±0 . 4 

13+0.4 

12+0.4 

4 . 320  nm 

film 

with  432 

molecules 

Li+  ion, 

1st 

shell 

106+14 

97+12 

87+17 

130+22 

Li ‘  ion , 

2nd 

shell 

10+0.3 

9+0.3 

10+0.3 

9+0.3 

X  ion. 

1st 

shell 

33±2 . 7 

16+1.2 

14+0.8 

11+0.6 

X  ion. 

2nd 

shell 

15+0.6 

16+0.6 

16+0.6 

15+0.5 
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1  able  III  Summary  ol  density  piolile  icsults  n|  halide  ions  in  two  water  films.  All 
positions  relative  to  closest  wall  at  \/\  U.‘>'4  !hu  film)  and  |/j  I.KOSnm 

(lor  4  ’Odum  liinti 


Li  +  a 

F 

Cl 

lit 

I  ~ 

2 . 362nm 

f  i  1m 

Znrar 

Zhb 

+0.084 

-0.084 

+0 . 034 

+0.058 

+0.058 

Z  pk 

Zhc 

+0.225 

-0.273 

-0.225 

-0.036 

-0.013 

^far 

Zhd 

+0 . 580 

-0.556 

-0 . 500 

-0.320 

-0.249 

4 . 32Qrtm 

film 

Zncar 

Zhb 

+0.044 

-0.007 

+0.062 

+0.062 

+0.062 

^pk 

Zhc 

+0.202 

-0.  176 

-0.044 

-0.017 

-0.017 

^far 

Zhd 

+0.414 

-0.440 

-0.3.35 

-0.308 

-0.202 

a  lithium  from  LiF  simulation 
b  /  .  point  of  nearest  approach  to  electrode 

c  /.  k,  position  of  peak  in  the  density 
d  /  position  of  farthest  point  from  electrode 
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Table  IV.  Summary  of  water  density  profile  results  for  1  iX  (X  =  I  ,(  I.Br.l )  two  water 
films.  All  positions  measured  in  nm  relative  to  nearest  hard  wall  at  |z„|  -=  0.934«ni  (for 
2.362nin  film)  and  \zh\  =  1.808  (for  4.320nm  film).  The  number  in  parentheses  is  the 
value  in  nm'1  of  the  density  profile 


LiF 

LiCl 

LiBr 

Lil 

2.362nm  film 

Peak  on 

-0.060 

-0.036 

-0.036 

-0.036 

anion  side 

(147) 

(143) 

(142) 

(135) 

Minimum  on 

-0.273 

-0.273 

-0.296 

-0.273 

anion  side 

(103) 

(  95) 

(  93) 

(  91) 

Peak  on 

0.036 

0.036 

0.036 

0.036 

cation  Side 

(133) 

(140) 

(138) 

(138) 

4.320nm  film 

Peak  on 

-0.044 

-0.044 

-0.044 

-0.044 

anion  side 

(157) 

(150) 

(152) 

(141) 

Minimum  on 

-0.255 

-0.282 

-0.255 

-0.255 

anion  side 

(  96) 

(  93) 

(  90) 

(  87) 

Peak  on 

0.044 

0.044 

0.017 

0.017 

cation  Side 

(146) 

( 149) 

(151) 

(147) 

FIGURE  CAPTIONS 


Figure  l.  Ion-vvatcr  pair  correlation  functions  for  lithium  and  the  halide  ions.  The 
four  correlation  functions  for  the  lithium  ions  are  superimposed  showing 
their  essential  equivalence.  For  chloride,  bromide,  and  iodide  the  second 
peak  near  0.55  nm  arc  also  essentially  the  same  within  numerical  error. 
The  dip  at  0.82  nm  indicated  by  the  arrows  was  due  to  the  truncation 
function. 

Figure  2.  Time  dependence  of  the  z  coordinate  of  the  lithium  ion  and  a  fluoride  ion 
during  a  simulation  lasting  1000  ps.  At  time  t  =  0  ps  both  ions  start  out 
near  the  cathode.  It  takes  the  anion  about  100  ps  to  migrate  across  the 
gap.  Note  that  the  electric  fields  keep  the  ions  separated  and  restricted  to 
vicinity  of  their  respective  electrodes.  The  ions  arc  not  free  to  diffuse 
throughout  the  gap.  Dashed  horizontal  lines  define  the  beginning  of  the 
hard  wall  repulsive  region  of  the  potential  starting  at  |z|  =  0.934  nm.  This 
is  effectively  the  boundary  surface  of  the  electrode. 

Figure  3.  Distribution  of  Li^  ,  F  ,  and  214  water  molecules  in  a  strong  electric  field 
(~2xl0?  V/cm)  across  the  gap  between  the  charged  plates.  Wall  potential 
shown  for  reference.  Circles  arc  Leonard-.! ones  radii  for  the  ions  and  water. 
In  this  simulation  the  anion  never  contact  adsorbs.  Contrast  this  with  the 
results  of  the  simulation  for  LiCl,  Lillr,  and  Lil.  Note  the  shallow  minimum 
in  the  water  distribution  near  the  average  position  of  the  anion. 

Figure  4.  Distribution  of  Li  * ,  F ,  and  430  water  molecules  in  a  strong  electric  field 
(~2xl07  V/cm)  across  the  gap  between  the  charged  plates.  Wall  potential 
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shown  for  reference,  Circles  arc  Lcnnard-Joncs  radii  for  the  ions  and  water. 
In  this  simulation  the  anion  never  contact  adsorbs.  Contrast  this  with  the 
results  of  the  simulation  for  LiCl,  Lilir,  and  Lil.  Note  the  shallow  minimum 
in  the  water  distribution  near  the  average  position  of  the  anion. 


Figure  5.  Time  dependence  of  the  z  coordinate  of  the  lithium  ion  and  a  a  chloride 
ion  during  a  simulation  lasting  1000  ps.  At  time  t  =  0  ps  both  ions  start 
out  near  the  cathode.  It  takes  the  anion  less  than  100  ps  to  migrate  across 
the  gap.  Note  that  the  electric  fields  keep  the  ions  separated  and  restricted 
to  vicinity  of  their  respective  electrodes.  The  ions  arc  not  free  to  dilfusc 
throughout  the  gap.  Dashed  horizontal  lines  define  the  hard  wall  region 
of  the  potential  starting  at  |z|  =  0.934  nm.  This  is  effectively  the  boundary 
surface  of  the  electrode.  Unlike  fluoride  the  chloride  ion  docs  touch  the 
boundary,  and  in  this  model  is  able  to  contact  adsorb  in  weak  fashion. 

Figure  6.  Distribution  of  LiH  ,  Cl,  and  214  water  molecules  in  a  strong  electric  field 
(~2xl07  V/cm)  across  the  gap  between  the  charged  plates.  Wall  potential 
shown  for  reference.  Circles  are  Lcnnard-Joncs  radii  for  the  ions  and  water. 
In  this  simulation  the  anion  makes  weak  contact  adsorption.  Contrast  this 
with  the  results  of  the  simulation  for  LiBr  and  Lil,  where  the  contact 
adsorption  is  more  frequent  and  therefore  stronger.  Note  the  shallow 
minimum  in  the  water  distribution  near  the  average  position  of  the  anion. 

Figure  7.  Distribution  of  Li  * ,  Cl  ,  and  430  water  molecules  in  a  strong  electric  field 
(~2xl07  V/cm)  across  the  gap  between  the  charged  plates.  Wall  potential 


shown  for  reference.  Circles  arc  Lcnnard-Joncs  radii  for  the  ions  and  water. 
In  this  simulation  the  anion  makes  weak  contact  adsorption.  Contrast  this 
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with  the  results  of  the  simulation  for  LiBr  ami  J.il,  where  the  contact 
adsorption  is  more  frequent  and  therefore  stronger.  Note  the  shallow 
minimum  in  the  water  distribution  near  the  average  position  of  the  anion. 


Figure  8.  Time  dependence  of  the  z  coordinate  of  the  lithium  ion  and  a  bromide  ion 
during  a  simulation  lasting  1000  ps.  Note  that  the  electric  fields  keep  the 
ions  separated  and  restricted  to  vicinity  of  their  respective  electrodes.  The 
ions  arc  not  free  to  diffuse  throughout  the  gap.  Dashed  horizontal  lines 
define  the  repulsive  region  of  the  wall  potential  starting  at  |z|  =  0.934  nm. 
This  is  effectively  the  boundary  surface  of  the  electrode.  Note  that  the 
bromide  ion  frequently  touches  the  boundary,  and  so  in  this  model  is  a 
moderate  contact  adsorber. 

Figure  9.  Distribution  of  Li +  ,  Dr  ,  and  214  water  molecules  in  a  strong  electric  field 
(~2xI07  V/cm)  across  the  gap  between  the  charged  plates.  V'all  potential 
shown  for  reference.  Circles  are  Lcnnard-Joncs  radii  for  the  ions  and  water. 
In  this  simulation  the  anion  engages  in  moderate  contact  adsorption. 
Contrast  this  with  the  results  of  the  simulation  for  LiF  where  there  is  no 
contact  and  Lil  where  the  iodide  makes  strong  contact.  Note  the  shallow 
minimum  in  the  water  distribution  near  the  average  position  of  the  anion. 

Figure  10.  Distribution  of  Li  +  ,  Br ,  and  430  water  molecules  in  a  strong  electric  field 

7 

(—2x10  V/cm)  across  the  gap  between  the  charged  plates.  Wall  potential 


shown  for  reference.  Circles  arc  Lcnnard-Joncs  radii  for  the  ions  and  water. 
In  this  simulation  the  anion  engages  in  moderate  contact  adsorption. 
Contrast  this  with  the  results  of  the  simulation  for  LiF  where  there  is  no 
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contact  and  Li  I  where  the  iodide  makes  strong  contact.  Note  the  shallow 
minimum  in  the  water  distribution  near  the  average  position  of  the  anion. 


Figure  11.  Time  dependence  of  the  z  coordinate  of  the  lithium  ion  and  a  iodide  ion 
during  a  simulation  lasting  400  ps.  Note  that  the  electric  fields  keep  the 
ions  separated  and  restricted  to  vicinity  of  their  respective  electrodes.  The 
ions  are  not  free  to  diffuse  throughout  the  gap.  Dashed  horizontal  lines 
define  the  hard  wall  region  of  the  potential  starting  at  |z|  =  0.934  nm. 
This  is  effectively  the  boundary  surface  of  the  electrode.  Note  that  the 
iodide  ion  spends  almost  all  its  time  at  the  boundary,  and  in  this  model  is 
classified  as  a  strong  contact  adsorber.  It  makes  only  short  time  excursions 
away  from  the  surface. 

Figure  12.  Distribution  of  Li+,  f,  and  214  water  molecules  in  a  strong  electric  field 
(~2xl07  V/cm)  across  the  gap  between  the  charged  plates.  Wall  potential 
shown  for  reference.  Circles  are  Lennard-Joncs  radii  for  the  ions  and  water. 
In  this  simulation  the  anion  exhibits  strong  contact  adsorption.  The  density 
profile  is  sharply  peaked  very  close  to  the  repulsive  wall.  Contrast  this  with 
the  results  of  the  simulation  for  LiF  where  there  is  no  contact.  Note  the 
shallow  minimum  in  the  water  distribution  near  the  average  position  of  the 
anion. 

Figure  13.  Distribution  of  Li  * ,  I  ,  and  432  water  molecules  in  a  strong  electric  field 

7 

(~2xl0  V/cm)  across  the  gap  between  the  charged  plates.  Simulation  time 


lOOOps.  Wail  potential  shown  for  reference.  Circles  are  Lennard-Joncs  radii 
for  the  ions  and  water.  In  this  simulation  the  anion  exhibits  strong  contact 
adsorption.  The  density  profile  is  sharply  peaked  very  close  to  the  repulsive 
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wall.  Contrast  this  with  the  results  ol  the  simulation  (or  Li!  where  there 
is  no  contact  Note  the  shallow  minimum  in  the  water  distribution  near 
the  average  position  of  the  anion. 
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